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OO . Abstract 

O 1 

t^j- ' LEPTO event generator is modified to describe the azimuthal modulations arising 

from Cahn and Sivers effects. The comparisons with some existing data in the 
current fragmentation region of SIDIS are presented. The predictions for Cahn and 
Sivers asymmetries in the target fragmentation region are presented for SIDIS of 12 
O ■ GeV electrons off proton target. 
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■ 1 Introduction 
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In P the role of parton intrinsic motion in semi-inclusive DIS (SIDIS) processes within 
QCD parton model has been considered at leading order; intrinsic is fully taken into 
account in quark distribution functions and in the elementary processes as well as the 
hadron transverse momentum, pj_, with respect to fragmenting quark momentum, see 
Fig.^. The average values of k± for quarks inside protons and p± for final hadrons inside 
the fragmenting quark jet where fixed by a comparison with data on Cahn effect [2] - the 
dependence of the unpolarized cross section on the azimuthal angle between the leptonic 
and the hadronic planes. The single spin asymmetry (SSA) A s ^ ^ recently observed 
by HERMES Collaboration j3] was successfully described by Sivers mechanism It 
was also shown that the Sivers distribution functions resulting from this analysis are 
compatible with the preliminary data from COMPASS collaboration 

The description of SIDIS within standard QCD parton model approach using the 
distribution and fragmentation functions is valid only in the current fragmentation region, 
CFR (xp > 0) and at high energies. A more general approach allowing to describe SIDIS 
in the whole kinematic region available for final hadrons is based on the LUND string 

1 In the following the notations of are used. 
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Figure 1: Three dimensional kinematics of the SIDIS process. 




fragmentation model and is incorporated into LEPTO event generator jHj. In the simplest 
case, corresponding to LO approximation of parton model, event generation in LEPTO 
proceeds in several steps: 

1. The active quark inside the nucleon is chosen according to the quark density function 
/ 9 (x,Q 2 ), 

2. The hard scattering kinematics is generated, 

3. The transverse momentum of the final quark is simulated with Gaussian k± and 
flat if distributions. Note that the transverse momentum of the final final quark is 
equal to that of initial quark for leading order hard subprocesses. 

4. The string fragmentation machinery of JETSET program [2j is applied to form the 
final hadrons. 

Within this approach the SIDIS cross section at LO can be expressed as 

dxdQ*dx F <PP T = ^ Hi /N (^F,k ± ,P T ), (1) 

where do ^n 2 ^ is the lepton-quark hard scattering cross section and H^ N (x, Xp, k±, P T ) 
represents the hadronization function of the system formed by struck quark with trans- 
verse momentum and target remnant. In the standard version of LEPTO the quark 
distribution function and the LO lepton-quark cross section are independent of ip, and, 
thus the final quarks are uniformly distributed in azimuthal angle. For final hadrons this 
implies also a uniform azimuthal distribution. However, already in unpolarized SIDIS the 
observed azimuthal distribution of hadrons is not flat. 
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In this paper two types of azimuthal modulation at the quark level and their influence 
on the produced hadron azimuthal distribution will be considered: 

• azimuthal modulation of the hard scattering cross section in unpolarized SIDIS 
(Cahn effect) 

• azimuthal modulation of the initial quark transverse momentum in SIDIS of unpo- 
larized leptons off the transversely polarized nucleon (Sivers effect). 

It is possible to incorporate both effects in the LEPTO event generator and obtain pre- 
dictions for azimuthal asymmetries in the whole kinematical region for the final hadrons. 
The way how the LEPTO code is modified to include Cahn and Sivers effects are described 
in Sec. |21and Sec. El respectively. In Sec. EJsome discussion and conclusions are presented. 



2 Including Cahn effect in LEPTO 

The Cahn effect j2j is a kinematical effect arising due to the presence of nonzero intrinsic 
transverse momentum of quarks in the nucleon. The underlying physics is very simple. 
The lepton-quark scattering cross section is given by the QED expression 

da e ^ ocs 2 + u 2 . (2) 

In the general case of non collinear kinematics Mandelstam variables depend on the quark 
transverse momentum and its azimuthal angle and at order 0(k±/Q) one has 

s 2 = 

u 2 = ^(1-^(1-4^^=1. (3) 
Then, the lepton-quark elastic scattering cross section is given by 




j~tq-+tq i ( 2 + y)V 1 -yk± ( s 

da xl ~ i + q 008 ^ (4) 

Eq. shows that the azimuthal angle of the final quark (and of the string's end asso- 
ciated with the struck quark) is now modulated with amplitude depending on y, Q and 
k ± . 

This effect can be introduced in the LEPTO event generator at the step 3) of the event 
generation, when the transverse momentum and azimuthal angle of the scattered quark 
are generated. To do this the generation of the quark transverse momentum, k±, is 
left unchanged and then the azimuthal angle is generated according to Eq. (@J). This 
leads to azimuthal modulation of the string axis (axis z on Fig. The momentum 
conservation means that the transverse momentum of the quark is balanced by that of 
the target remnant, which in turn means that the azimuthal angle of the target remnant 
fqq = (p + 7i. Then, one expects that the azimuthal angle of the hadrons in the target 
fragmentation region (TFR), xf < 0, will be modulated with shifted a phase by n with 
respect to that in CFR. 
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Figure 2: The xp dependence of (cos <f>h) jw\ (y) for charged hadrons compared with EMC data. 







Figure 3: Predictions of modified LEPTO for xp dependence of (cos^/j) for different hadrons 
produced in 12 GeV unpolarized SIDIS process. 



Data on azimuthal dependencies of SIDIS covering a large xp range have been ob- 
tained by the EMC Collaboration [S] for a beam energy of 280 GeV. The xp depen- 
dence of (cos0/j) /wi(y), where Wi(y) = (2 — y)y/l — y/ (1 + (1 — y) 2 ), obtained by using 
modified LEPTO for EMC kinematics are presented in Fig. |2] together with data points 
from |S]. The simulations has been done with LO setting of LEPTO (LST(8)=0) and with 
values of the parameters describing intrinsic kp (PARL(3)=0.5) and fragmentation pp 
(PARL(21)=0.45) as adopted in 0. 
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The predictions of modified LEPTO for (cos^) of different hadron (ir + , 7r~, 7r° and 
p) produced in SIDIS on a proton target at future CEBAF 12 GeV facility at JLab [Hj 
are presented in Fig. |3| One can see from Fig. |2] and Fig. |3] that the predicted mean 
value of cos 0^ in the CFR is negative (cos (ph)cFR < 0, while in the TFR is positive 
(cos 4>h)TFR > 0, as suggested by arguments based on transverse momentum conservation. 

3 Including Sivers effect in LEPTO 

The azimuthal modulation of the string transverse momentum in the previous section 
was due to Cahn effect - the dependence of the non planar lepton-quark scattering cross 
section on the quark azimuth. The quark distribution, f q (x,k±) itself is independent of 
quark azimuthal angle. 

The situation is different when one considers SIDIS on a transversely polarized nu- 
cleon. Now a correlation between transverse momentum of quark and target transverse 
polarization is possible - the so called Sivers effect jlj. For quite some time it was be- 
lieved that this correlation is forbidden because of T-invariance of the strong interactions. 
However the spectator model calculations ^U] demonstrated that there exists a nonzero 
SSA in SIDIS when the final state interaction between struck quark and target remnant 
is taken into account. Then, the effective description of this SSA is possible within QCD 
factorized approach by introducing a new distribution function - the Sivers function (for 
further discussion see, for example, [TT]). 

The unpolarized quark (and gluon) distributions inside a transversely polarized proton 
(generically denoted by p\ with p^ denoting the opposite polarization state) can be written 
as: 

/,/pt (x, k±) = f q / p (x, k ± ) + - A N f q/p i (x, fcj.) S T ■ (P x k±) , (5) 

where P and S T are respectively the proton momentum and transverse polarization vec- 
tor, and k± is the parton transverse momentum; transverse refers to the proton direction. 
Eq. dHJ) implies 

f q / p t(x, k ± ) + f q /pi(x, fcjj = 2f q / p (x, k±) , 

/,/pT (x, fex) - fq/pi (x, k±) = A N f q/p i (x, fc x ) S T ■ (P x fe ± ) , (6) 

where f q / p (x,k±) is the unpolarized parton density and A N f q / p i (x, k±) is referred to as 
the Sivers function. Notice that, as requested by parity invariance, the scalar quan- 
tity S T ■ (P x k±) singles out the polarization component perpendicular to the P — k± 
plane. For a proton moving along — z and a generic transverse polarization vector 
S T = \S T \ (cos 4>s, sin <p s , 0) (see Fig. [TJ) one has: 

S T ■ (P x fcj.) = \S T \ sin(v9 - (p s ) = \S T \ sin0 5 ^ , (7) 

where (<p — <fis) = 4>siv is the Sivers angle. 

In the Sivers function for each light quark flavor q = u,d are parameterized in the 
following factorized form: 

A N f q/p r(x, k ± ) = 2N q {x) h(k ± ) f q/p (x, k ± ) , (8) 
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where 



N q {x) = JV fl s- (1 - x) b « , (9) 

a q b q 

HkJ^VTe^e-^, (10) 

where N q , a q , b q and M (GeV/c) are parameters. Then Eq. ©can be rewritten as 

f q / p i(x,k±) = f q / p [x,k±)(l + \S T \J\f q (x)h(k±) sm4> S iv\- (11) 

Again, the Sivers effect is incorporated into LEPTO at the stage 3) of the event gener- 
ation in the same way as for the Cahn effect but now the azimuthal angle is generated 
according to Eq. (flljl . For simulations the following set of parameters compatible with 
those obtained in [T] have been used: N u = = 0.5, Nd = = —0.2, a q = 0.3, b q = 2 
and M 2 = 0.36 (GeV/c) 2 . 

In Fig.|3]the results of simulation for HERMES experimental conditions are compared 
with observed Sivers asymmetries [3]. 
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Figure 4: HERMES data on A S J^^ jSj for scattering off a transversely polarized proton 
target. The curves are the results of simulations obtained with modified LEPTO. 



Future facilities as Electron Ion Colliders or upgraded JLab will have larger kinematic 
coverage and will offer the possibility of studying the Sivers effect also with hadrons 
produced in the TFR. As an example, the simulations have been done for 12 GeV electron 
SIDIS of a proton target. The DIS cut Q 2 > l(Gev/c) 2 and W 2 > AGev 2 and a cut on the 
produced hadron transverse momentum Py > 0.05 GeV/c was imposed. The predictions 
for Xp, x and Pt dependencies for JLab kinematics are presented on the Fig. Fig. El 
and Fig. [7[ respsctively. 
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Figure 5: Predicted dependence of At™^ Ph on xj? for different hadrons produced in SIDIS 
of 12 GeV electrons off a transversely polarized proton target. 
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Figure 6: Predicted dependence of A V j, 



sin(tp h -(p s ) 



on x for different hadrons produced in the TFR 



(xf < —0.1) of SIDIS of 12 GeV electrons off a transversely polarized proton target. 
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Figure 7: Predicted dependence of Aj^f Vh on for different hadrons produced in the 
TFR (xp < —0.1) of SIDIS of 12 GeV electrons off a transversely polarized proton target. 

4 Discussion and Conclusions 

In this article the way of modifying the standard LEPTO event generator in order to 
include the azimuthal asymmetries arising from Cahn and Sivers effects is described. 
Only LO effects have been taken into account. The azimuthal modulations for Cahn 
and Sivers effects have different origins. In the case of Cahn effect the initial quark 
transverse momentum is independent of azimuthal angle but the hard scattering cross 
section in a non planar kinematics depends on the final quark azimuthal angle. In the 
case of Sivers effect already the initial quark transverse momentum has an azimuthal 
modulation. The azimuthal asymmetries are introduced in both cases by changing the 
struck quark/string azimuthal distribution during event generation. The hadronization 
part of program (JETSET) is left unchanged. The possible influence of the higher twist 
distribution functions as well as possible modifications of hadronization in the case of 
polarized target ^2] have been ignored. 

The advantage of this MC based approach compared to standard QCD factorized 
approach is in the full coverage of produced hadron phase space. The modified generator 
will be useful for complete MC simulations of experiments including Cahn and Sivers 
effects both in the CFR and in the TFR and also for global analysis of these effects. 

Figs. 121 and 0] demonstrate that the modified LEPTO event generator well describing 
the data in the CFR both for Cahn and Sivers asymmetries. The description of Cahn 
effect in the TFR looks unsatisfactory. This discrepancy can be explained either by some 
unaccounted contributions in the TFR or by insufficient precision of experimental data. 
One can notice from the experimental points in Fig.|2]that the integrated value of (cos^) 
for charged hadrons in the CFR is not compensated by that in TFR. It seems improbable 
that this imbalance can be compensated by larger values of (cos0/ 4 ) of neutral hadrons 
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at xf — — 1. 

The new high statistic measurements will allow to check the predictions of the ap- 
proach presented here and better understand the effects of the quark intrinsic transverse 
momentum and hadronization mechanism in SIDIS. 
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